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ABSTRACT 

Recent timing observations of PSR J0045-7319 reveal that the neutron star/B 
star binary orbit is decaying on a time scale of |-P rb/-forb| = 0.5 Myr, shorter than 
the characteristic age (r c = 3 Myr) of the pulsar (Kaspi et al. 1996a). We study 
mechanisms for the orbital decay. The standard weak friction theory based on static 
tide requires far too short a viscous time to explain the observed P rb- We show 
that dynamical tidal excitation of g-modes in the B star can be responsible for the 
orbital decay. However, to explain the observed short decay timescale, the B star must 
have some significant retrograde rotation with respect to the orbit — The retrograde 
rotation brings lower-order g-modes, which couple much more strongly to the tidal 
potential, into closer "resonances" with the orbital motion, thus significantly enhancing 
the dynamical tide. A much less likely possibility is that the g-mode damping time is 
much shorter than the ordinary radiative damping time. The observed orbital decay 
timescale combined with a generic orbital evolution model based on dynamical tide can 
be used as a "timer", giving an upper limit of 1.4 Myr for the age of the binary system 
since the neutron star formation. Thus the characteristic age of the pulsar is not a 
good age indicator. Assuming standard magnetic dipole braking for the pulsar and no 
significant magnetic field decay on a timescale ^ 1 Myr, the upper limit for the age 
implies that the initial spin of the neutron star at birth was close to its current value. 



Subject headings: binaries: close - pulsars: individual (PSR J0045-7319) - stars: 
neutron - stars: oscillations - stars: rotation - hydrodynamics - supernovae 



1. Introduction 

One of the fundamental questions in the study of pulsars concerns the initial conditions of 
neutron stars at birth. In particular, the initial spin of neutron star is related to such issues 
as angular momentum coupling in the progenitor stars, supernova explosion mechanisms and 
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gravitational wave generation from core collapse and nascent neutron stars. Unfortunately, except 
for the Crab pulsar, for which the initial spin can be directly inferred from the measured pulsar 
period P p and its time-derivatives P p , P p and the known age, our knowledge about this quantity 
is rather limited. Statistical studies of the evolution of pulsar population (e.g., Narayan 1987), the 
energetics of pulsar nebulae (Helfand & Becker 1987) and possible old pulsar /supernova remnant 
associations (Kaspi et al. 1996b) give some indications that neutron stars are formed rotating at 
a moderate spin period 10 — 100 ms. However, all these methods suffer from uncertainties, and 
other independent constraints on the initial spins of neutron stars are highly valuable. 

The PSR J0045-7319 binary contains a radio pulsar (P p = 0.93 s) and massive B-star 



companion in an eccentric, 51.17 days orbit ( Kaspi et al 1994] ). This system has recently been 



shown to exhibit spin-orbit precessions due to the rapid, misaligned rotation of the B star, which 



strongly suggests that the neutron star received a kick at birth from asymmetric supernova (Lai 



Bildsten & Kaspi 1995, Kaspi et al. 1996a). Recent timing observations also reveal that the orbit 
is decaying on a time scale |-P rb/-forb| = 0.5 Myr, shorter than the pulsar's characteristic age 
t c = Pp/i^Pp) = 3 Myr ( |Kaspi et al. 1996a] ) . In this paper, we study the physical mechanisms for 
the rapid orbital decay, and demonstrate the potential of using the orbital decay time to constrain 
the age and initial spin of the pulsar. 

The fiducial numbers we adopt for the current PSR J0045-7319 system are: pulsar mass 
M p = IAMq, companion mass M c = 8.8M@, radius R c = 6ARq, and orbital period P or b = 51.17 
days, semi-major axis a = 20R C , eccentricity e = 0.808 ( Kaspi et al. 1996a , Bell et al 1995|) . 



2. Mechanisms for the Orbital Decay 



Since the mass loss from the B star is completely negligible ( |Kaspi et al. 1996a ), tidal 



interaction is the most natural cause for the orbital decay in the PSR J0045-7319 system. 
There are two kinds of tides: static and dynamical. The static tide corresponds to the global, 
quadrupolar distortion of the star. In the presence of finite viscosity, a tidal lag angle #i ag oc t~^ c 
develops between the tidal bulge and the line joining the pulsar and the B star. The orbital 
evolution and spin evolution are driven by the viscous tidal torque ~ {GM p R^ c / cfi)9\. Ag . However, 
this "weak friction theory" based on static tide cannot explain the orbital decay in the PSR 
J0045-7319 binary for two reasons: (i) Using the general tidal equations (e.g., Alexander 1973) 
together with constraints on the stellar rotation rate and spin-orbit inclination angle, we have 
found that a viscous time t v i sc of less than 30 days is needed in order to explain the observed -P or b; 
Such a small viscous time is almost certainly impossible, (e.g., even if the star were completely 
convective, the viscous time would still be longer than one year), (ii) The same amount of torque 
that induces the orbital decay also acts on the B star; Since the stellar spin S is much less than 
the orbital angular momentum L, the static tide would lead to rapid spin-orbit synchronization 
and alignment if it were to explain the observed P rb, m contradiction with observations. In the 
following, we discuss the physics and outline the main result of a dynamical tide theory that can 
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explain the observations. A more complete analyses will be presented in a future paper. 

At each periastron passage, the pulsar excites non-radial (mainly I = 2 quadrupolar) 
oscillations in the B star, therefore energy and angular momentum are transferred between the 
orbit and the star. Write the Lagrangian displacement in the rotating frame of the star as 
£(r, t) = J2a a a(t)£,a(r), where £ Q (r) is the eigenmode, normalized via Jd 3 xp£* ■ £ a = 1, and 
a = {n\m} is the mode index (n is the radial order of the mode, A and m specify the angular 
eigenf unction; note that A = I only for nonrotating star). The mode amplitude a a (t) is then 
governed by the equation 

d Q + Jtfla + 2j a d a = (1) 

where iv a is the mode angular frequency (in the rotating frame), 7 Q the amplitude damping rate, 
D(t) is the orbital separation, the orbital phase, Q s is the rotation rate of the companionP], 
Wi m is a constant as defined in Press and Teukolsky (1977), and Q a is the tidal coupling coefficient 
as given by Q a = fd 3 xp£l ■ X7(r l Yi m ). 

The energy transfer between the star and orbit during each periastron passage depends on the 
phases of the oscillation modes, thus in general varies from one passage to another (e.g., Kochanek 
1992; Mardling 1995). An important measure of the strength of dynamical tide is the energy 
AE S transferred to the star (and the corresponding angular momentum transfer A J s ) during the 
"first" periastron when there is no oscillation present initially. This can be calculated using the 
formalism of Press and Teukolsky (1977), generalized to include stellar rotation and elliptic orbit. 
To the leading (1 = 2 quadrupole) order, we have 

GM? /P.\* 

where rj = (M c /M t ) 1 ^ 2 (a p /R c ) 3 ^ 2 is the ratio of the time for periastron passage and the stellar 
dynamical time, Mt = M p + M c is the total mass, and a p = a(l — e) is the periastron distance. 
The factor F a = 1 — (i/u a ) J d?x • (O s x £ Q ) is of order unity, and K[ m ^ depends on the orbital 
trajectory and the mode frequency: 

W lm fPorW^ 1 
-^orb/2 

with the integration centered around periastron. We can obtain a similar expression for the 
angular momentum transfer: 




Wi m f p ° A / 2 1 

Klm,a = TT~ / rf *7v+T COS + ~ mQ, s t], (3) 

27T J-P„ rh /2 D L + i 



1 The observation of orbital precession implies that the stellar spin and the orbital angular momentum are not 
aligned. Here we assume they are aligned for simplicity. Thus Q a should be understood as the component of the spin 
rate perpendicular to the orbital plane. 
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Clearly, the most strongly excited modes which contribute most to the energy transfer are those 
(i) propagating in the same direction as the orbit (the m = —2 modes); (ii) having frequencies in 
the inertial frame (u> a — mQ s ) ~ (— m<3? p ), where <& p = f2 orb (l + e) 1 / 2 /(l — e) 3 / 2 = 27r/(3.2 days) is 
the orbital frequency at periastron (" resonance" condition); and (hi) having relatively large Qc*. 
For nonrotating star, these modes are g-modes of radial order n = 5 — 9 (i.e., g§ — gg). 

The rapid rotation of the B star changes its g-modes significantly. Since the spin frequency 
is comparable to the relevant g-mode frequencies, a perturbation treatment of the rotation effect 
is not valid. Neglecting the centrifugal force, the g-modes of rotating star can be calculated using 
the "traditional approximation" (e.g., Unno et al. 1989), where the Coriolis force associated with 
the horizontal component of the spin angular velocity is ignored (valid because the horizontal 
displacement of the g-mode is larger than the radial displacement). We correct the g-mode 
frequencies of a T = 4/3, T\ = 5/3 nonrotating polytropef] using the numerical result of 
Berthomieu et al. (1978) to obtain the mode frequencies of the rotating star. Since |w a /fJ s | is 
greater than unity (although not much greater) for the modes that are strongly excited, the 
rotational corrections to the mode frequencies and wavefunctions are not large (less than 40% in 
uj a ). We therefore adopt Q a to be approximately the same as that of the non-rotating star; its 
validity is also indicated by an asymptotic analysis (Rocca 1987). 

Figure 1 shows the function T2(r?) for different values of £l s . We see that rotation can 
significantly change the strength of dynamical tide. In particular, at n = 7 (appropriate for 
the current PSR J0045-7319 system), retrograde rotation (f2 s < 0) can increase AE S from the 
nonrotating value by more than two orders of magnitude. Physically, this dramatic effect of 
rotation on the tidal strength comes about because (i) the retrograde rotation brings lower-order 
g-modes, which couple to the tidal potential more strongly, into closer "resonances" with the 
orbital motion; (ii) The tidal coupling coefficient Q a depends very strongly on the order of the 
mode. For example, for Cl s = O s (GM c /^) -1 / 2 = —0.4, the most strongly excited modes are 
94,-95, while g\ — 53 are also much more strongly excited compared to the nonrotating case. The 
function £2(7?) has similar behavior as 22(77), and we find that typically S2 — 2T2 to within 10% 
(for \(l s \ < 0.5). 

The radiative damping times of f-mode and low-order g-modes of a massive main sequence 
star range from 10's to 100's of years (Saio & Cox 1980; Dziembowski et al. 1993), much longer 
than -P rb- Figure 2 shows an example of the dynamical tidal evolution obtained by integrating 
Eq. (1) and the orbital equations. The average mode energy settles into a constant value E s ~ AE S 
(a "fixed point" of the dynamical system) after a few damping times, and the orbital energy 
decreases secularly according to 

E orb ~ -2 7 |A£ orb [ = -2 7 (A£ S + ft s AJ s ), (5) 



2 This is appropriate for a massive star, although the g-modes depend somewhat on the size of the convective core, 
and hence on the stellar age. We have tested that |A_E s j is rather insensitive to the core size. 
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where 7 a = 7 = 1/idamp (i-e., we assume that the damping times are the same for different 
modes; this is approximately valid since only a small number of modes are strongly excited). With 
S2 — 2T2, we have AJ S ~ 2(R^/GM C ) 1 ^ 2 AE S , and the secular orbital decay rate is then given by 

Pprb i + 2A 5 / 100P orb \ / r 2 \ 
P orb ~ 0.24 Myr ^ i damp J Uo"V ' 

For nonrotating star, we find T2 ~ 2.6 x 1CP 4 , thus a very short damping time tdamp — 5.5P rb is 
needed to explain the observed P rb- However, when the star has significant retrograde rotation, 
the observed orbital decay rate can be explained with ordinary radiative damping. 

If the mode damping time is shorter than the orbital period, then the secular orbital 
decay rate should be given by P Q rb ~ — | A£7 or b | / [ £?orb I • This would require significant 
enhancement to the mode dissipation compared to ordinary radiative damping. Nonlinear 
mode damping (Kumar & Goodman 1996) is unimportant since the mode energy 
E S /(GM 2 I R c ) ~ 10 _7 (tdamp/100P O rb) ) implied from the observed orbital decay rate, is 
much smaller than the threshold E s jh/(GM 2 /R c ) ~ (u) a tth)~ 0A ~ 1CP 3 needed for the three-mode 
parametric resonant coupling (where tth ~ 10 4 yr is the thermal time of the star, and we have 
approximated the damping time of the n ~ I >> 1 mode by tth /I 8 )- We consider this possibility 
rather unlikely. 



3. Orbital Decay Model and Constraints on the Pulsar Age and Initial Spin 

We now consider the long-term evolution of the binary system. Equations (2) and (5) 
provide a scaling relation for the orbital decay rate. Let ^2(77) oc r]~ 4u (with v ~ 1, 0.5, 0.2 for 
Cl s = 0, —0.2, —0.4 respectively), we obtain 

Por b = -AP~ 7 h /3 ~ 4U (l-e)-^\ (7) 

where the constant A can be fixed by the observed current P Q rb value. Using 
L = —j s ~ -2(Rl/GM c )^ 2 E s , we have 

L _ 1 P orb ,o(l - eg) 1 / 2 2P orb 

L (5 P orb (l - e 2)V2 3 P orb ' U 

where the subscript "0" indicates current values, and the parameter (3 ~ 50(1 + 2Q S ) 
ranges from 10 to 50 for —0.4 < Q s < 0. Note that the long-term orbital 
angular momentum change^] is much less than the change in the orbital energy since 



3 The orbital angular momentum is transferred into the spin, whose evolution can be easily evaluated: Since 
tl s /Q s ~ (0.04/f2 3 )(P or b/P O rb) (with the moment of inertia 0.1M C R^.), the timescale for changing the stellar spin 
(both its magnitude and direction) is longer than the orbital decay timescale. This feature is qualitatively different 
from static tide. 
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AJ S /L\ ~ |AE s / J E; orb |^ orb (GM c /i?3)- 1 /2( 1 _ e 2)-i/2 |A£ s /£ orb j. The evolution of 



parameters v and /? have been chosen to reflect their allowed ranges. The results are shown 
in Figure 3. We see that e — > 1 and P or b — > oo as (— t) approaches a fixed value <^ 1.4 Myr, 
independent of the parameters v and (3. This feature can be understood as follows: Since the 
angular momentum transfer in a periastron passage is relatively small compared to the energy 
transfer, the periastron distance a p and hence AE S remain approximately constant during the 
evolution. Thus |-P or b/-Porb| °c AE s /\E or t,\ is larger at earlier times. 

The observed P or b and our generic orbital evolution model therefore set an upper limit ~ 1.4 
Myr to the age of the binary system since the last supernova []. This is significantly smaller than 
the characteristic age r c of the pulsar, implying that the latter is not a good age indicator. In 
principle, the discrepancy between the age and r c may be explained in several ways: (i) The pulsar 
has been slowed down by an earlier phase of mass accretion; (ii) The pulsar braking index is much 
larger than 3 (the canonical dipole value); (iii) The pulsar magnetic field decays on a timescale 
<^ 1 Myr; (iv) The initial spin period of the pulsar is not much shorter than the current value. The 
extremely small mass loss from the B star (Kaspi et al. 1996a) argues strongly against (i). The 
observed braking indices measured for three of the youngest pulsars range from 2.0 to 2.8, and 
in view of the statistical evidence and physical argument against rapid field decay, we consider 
(ii) and (iii) to be rather unlikely. Using the standard magnetic dipole braking for the pulsar 
spin-down, the constraint on the pulsar age implies that the initial spin period of the pulsar is 
longer than 0.7 second, very close to its current value. Thus PSR J0045-7319 was formed rotating 
very slowly, and its progenitor must also have rather slow rotation. We note that such a small 
initial spin rate of the pulsar may result from the random torque on the nascent neutron star due 
to asymmetry in the supernova (e.g., Burrows et al. 1995). 



4 We assume there is no major structure change in the B star during the orbital evolution. This is most likely to 
be valid because the Kelvin-Holmholtz timescale (~ 10 4 years) is much shorter than the orbital evolution time. 

5 In the unlikely case of idamp Ss -Porb, the scaling relation for P or b is modified from Eq. (7) to P or b oc 
-• P orb° /3 ~ 4 "( 1 " e)- 6 < 1+I/) , and an absolute upper limit to the age cannot be obtained (see the dotted line in Fig. 3). 
However, the age can still be constrained since the probability of forming a highly eccentric, bound orbit after 
supernova is small. 



eccentricity is given by 




(9) 



For the current PSR J0045-7319 binary, we find e ~ 10 s , consistent with the observational 
upper limit for e (Kaspi et al. 1996a). 

Equations (7)-(9) can be integrated backward to give _P rb and e at earlier timesQ. The 
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4. Discussion 

Our analyses of PSR J0045-7319 binary orbital decay show that the B star companion is 
most likely to have a retrograde rotation with respect to the orbit. Since the spin of the B star 
is expected to have been aligned with the orbital angular momentum before the supernova, the 
current misaligned configuration could have come about only if the neutron star received a kick at 
birth due to asymmetry in the supernova. The kick velocity must (i) have had a component out 
of the original orbital plane in order to misalign L and S; (ii) have had a significant component 
in the direction opposite to the orbital velocity of the progenitor in order to reverse L. The 
current timing data yield two degenerate solutions for the range of the spin-orbit inclination angle, 
allowing for both prograde and retrograde rotation (Kaspi et al. 1996a). Long-term (i.e., some 
fraction of the precession period ~ 500 years) timing observation should distinguish these two 
possibilities. Dedicated optical observation of the companion may also give useful constraints on 
the excited modes (Kumar et al. 1995). 

Our analyses have also demonstrated that the orbital decay can be used to put concrete 
constraints on the pulsar age and initial spin. The tide-induced orbital decay of the PSR 
B1259-63/Be star binary (the only other known radio pulsar/main sequence star binary) is too 
slow to be observable owing to its larger orbital separation at periastron. Finding more systems 
similar to PSR J0045-7319 will allow for determination of systematic constraints on the physical 
conditions of neutron stars at birth and supernova characteristics. 

I thank Vicky Kaspi for sharing the observational data with me and helpful discussion and 
comments. I also thank Lars Bildsten, Peter Goldreich, Norm Murray and Yanqin Wu for useful 
discussion. This research is supported by the Richard C. Tolman Fellowship at Caltech, NASA 
Grant NAG 5-2756, and NSF Grant AST-9417371. 
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Fig. 1. — The function T2(rj) (denned in Eq. [2]) for different rotation rates: Cl s = (solid lines), 0.2 
(dotted lines), —0.2 (short-dashed lines), and —0.4 (long-dashed lines). Negative £l s corresponds 
to retrograde rotation. The heavy lines are for elliptic orbit with e = 0.808 (the current value for 
the PSR J0045-7319 system), and the light lines for parabolic orbit. 

Fig. 2. — Evolution of 5E or b = E or b — £? or £,,o and E s (mode energy in the rotating frame) due to 
dynamical tide (with ti, s = —0.2). Only the g% — gs modes are included in the calculation, and 
relatively large damping rates are chosen for clearer illustration: 7 = 0.1/P O rb (solid curves) and 
7 = 0.05/P or b (dotted curves). 

Fig. 3. — The long-term (backward) evolution of the binary period and eccentricity (t = —3 
Myr corresponds to the pulsar's characteristic age). The solid lines are for (3 = 10, with 
v = 0.2, 0.5, 1 (left to right), the short-dashed lines are for (v, (5) = (0.5,20), the long-dashed lines 
for (u, (3) = (0.5, 50). The dotted curves are for a model assuming tdamp ^ -Porb with (u, (3) = (1, 50) 
(see footnote 5). 
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